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A search for supersymmetry in the context of general gauge-mediated (GGM) 
breaking with the lightest neutralino as the next-to-lightest supersymmetric particle 
and the gravitino as the lightest is presented. The data sample corresponds to an 
integrated luminosity of 36 pb _1 recorded by the CMS experiment at the LHC. The 
search is performed using events containing two or more isolated photons, at least 
one hadronic jet, and significant missing transverse energy. No excess of events 
at high missing transverse energy is observed. Upper limits on the signal cross 
section for GGM supersymmetry between 0.3 and 1.1 pb at the 95% confidence level 
are determined for a range of squark, gluino, and neutralino masses, excluding 
supersymmetry parameter space that was inaccessible to previous experiments. 
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Supersymmetry (SUSY), in particular the version based on gauge-mediated breaking HHZl/ is 
of particular theoretical interest for physics beyond the standard model (SM). Supersymmetry 
stabilizes the mass of the SM Higgs boson, drives the grand unification of forces, and incorpo- 
rates dark matter candidates within its framework. Previous searches for SUSY with gauge- 
mediated breaking |8HT4| were performed using a minimal model [15 1 as a benchmark. In that 
model, in order to reduce the number of free parameters, several assumptions, including gaug- 
ino mass unification, are made. These assumptions lead to a mass hierarchy in which strongly 
interacting SUSY partners are much heavier than the lightest chargino and neutralino. For ex- 
ample, the current best lower limit on the neutralino mass [9] of 175 GeV corresponds to gluino 
and squark mass limits of well above 1 TeV. In the more general case, the masses of strongly 
interacting SUSY partners can be much lighter, leading to large production cross sections at the 
LHC, and allowing, even at low integrated luminosity, for the exploration of parameter space 
inaccessible at previous colliders. 

In this Letter we consider a general gauge-mediation (GGM) SUSY scenario I[i6l \Y7] . with the 
gravitino as the lightest SUSY particle (LSP) and the lightest neutralino as the next-to-lightest 
(NLSP). In the following, it is assumed that the neutralino decays promptly to a gravitino and a 
photon. Cases with either a large neutralino lifetime or a large branching fraction into a Z and a 
gravitino are not considered in this Letter. The gravitino escapes detection, leading to missing 
transverse energy (E™ ss ). If R-parity [18J is conserved, strongly-interacting SUSY particles are 
pair-produced. Their decay chain includes one or several quarks /gluons, and a neutralino, 
which in turn decays to a photon and a gravitino. The topology of interest for this search is, 
therefore, two or more isolated photons with large transverse energy (Ej), at least one hadronic 
jet, and large missing transverse energy E™ 1SS . 

A detailed description of the CMS detector can be found elsewhere fl9| . The detector's central 
feature is a superconducting solenoid providing a 3.8 T axial magnetic field along the beam di- 
rection. Charged particle trajectories are measured by a silicon pixel and strip tracker system, 
coveringO < <p < 2 n in azimuth and \r/\ < 2.5, where the pseudorapidity rj = — lntan0/2, and 
9 is the polar angle with respect to the counterclockwise beam direction. A lead-tungstate crys- 
tal electromagnetic calorimeter (ECAL) and a brass/ scintillator hadron calorimeter (HCAL) 
surround the tracker volume. For the barrel calorimeter ( 1 1] \ < 1 .479), the modules are arranged 
in projective towers. Muons are measured in gas detectors embedded in the steel return yoke of 
the magnet. The detector is nearly hermetic, allowing for reliable measurement of E™ 1SS . In the 
2010 collision data, photons with energy greater than 20 GeV are measured within the barrel 
ECAL with a resolution of better than 1%, which is dominated by intercalibration precision. 

The data used in this analysis, corresponding to an integrated luminosity of 36 pb _1 , were 
recorded during the 2010 run at a center-of-mass energy (\/s) of 7 TeV at the Large Hadron 
Collider (LHC) at CERN. A two-level trigger system was used, based on the presence of at 
least one photon with a minimum transverse energy of 30 GeV. This data sample is used for 
the selection of both signal candidates and control samples used for background estimation. 
The efficiency for offline selected events to pass the trigger is estimated to be above 99%. 

The photon candidates are reconstructed from clusters of energy in the ECAL. Candidates are 
required to have Ej > 30 GeV and \rj\ < 1.4. The ECAL cluster shape is required to be consis- 
tent with that expected from a photon, and the energy detected in the HCAL behind the photon 
shower is required not to exceed 5% of the ECAL energy. To suppress photons originating from 
quark/ gluon hadronization, the photons are required to be isolated from other activity in the 
tracker, ECAL and HCAL. The scalar sums of transverse energies of tracks and calorimeter de- 
posits within AR = y/ (A^) 2 + (A</>) 2 = 0.4 of the candidate's direction are determined, after 
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excluding the contribution from the candidate itself. These isolation sums are required to be 
< 0.001 x £ T + 2.0 GeV, 0.006 x £ T + 4.2 GeV, and 0.0025 x E T + 2.2 GeV, with E T in GeV, for 
the tracker, ECAL, and HCAL, respectively. 

Photons that fail either the shower shape or track isolation requirement are referred to as fake 
photons. Most of these fake photons are quarks/ gluons that have hadronized predominantly to 
7T°s decaying into photons. They are used for the background estimation derived from data. 

The criteria above are efficient for selection of both electrons and photons. To separate them 
reliably we search for hit patterns in the pixel detector consistent with a track from an electron 
(pixel match). The candidates with (without) pixel match are considered to be electrons (pho- 
tons) [20 1 . Events in the signal sample (referred to below as the 77 sample) are required to have 
at least two photon candidates. 

Jets are reconstructed from energy deposits in the calorimeters using the anti-fcj- clustering 
algorithm |2"T] with a size parameter of 0.5. The jet energy is corrected using reconstructed 
tracks p2| . Selected jets must have Ej > 30 GeV and \tj\ < 2.6. In addition, to separate real 
jets from anomalous HCAL signals, the fraction of energy contributed to the jet shower by the 
highest energy HCAL channels must be < 98%, the jet must have no single HCAL channel 
containing more than 90% of its total energy, and finally the ECAL energy fraction of the jet 
must be > 1%. To be retained in the signal sample, events must contain at least one such jet 
isolated from both of the two highest- Ej photon candidates by AR > 0.9. 

The E™ 1SS is determined using calorimeter energy deposits. Corrections are applied by replac- 
ing calorimeter tower energies matched to charged hadrons and muons with their correspond- 
ing charged-track momenta [23 1. 

Events are generated in the benchmark SUSY model Ij24"t using PYTHIA 6.4 II25II in a three- 
dimensional grid of the NLSP, gluino, and squark masses. The soft masses of the squarks 
are taken to be degenerate. Sleptons and all other gauginos except the NLSP are assigned 
a mass of 1.5 TeV. The QCD production cross section at NLO is calculated for these points 
using PROSPINO 2.1 11261 , and is dominated by gluino-gluino, gluino-squark, and squark-squark 
production. The generated events are passed through the CMS detector simulation program, 
which is based on the GEANT4 |27) package, and reconstructed using the same program as for 
the collision data to ensure that all features of the data, such as trigger and reconstruction, are 
applied to the Monte Carlo (MC) simulated SUSY signal sample. 

The SUSY signal can be mimicked in several ways. Irreducible backgrounds from SM processes 
such as Z(— > 1/17)77 and W(— > £v)jy are negligible. The main backgrounds arise from SM 
processes with misidentified photons and /or mismeasured E™ 1SS . The dominant contribution 
comes from mismeasurement of E™ ss in QCD processes such as direct diphoton, photon plus 
jets, and multijet production, with jets mimicking photons in the latter two cases. This back- 
ground is referred to as the QCD background. The strategy for determining this background is 
to use control samples that are kinematically similar to the candidate sample and that can be 
reasonably assumed to have no significant genuine E™ 1SS . Two such samples are used. The first 
is a sample containing two fake photons, referred to as the ff sample below, comprising QCD 
multijet events. The second (ee) sample contains events with two electrons |28| with invariant 
mass between 70 and 110 GeV, and is dominated by Z — > ee decays. 

The Ej resolution for electrons and fake photons is similar to the resolution for photons and is 
much better than the resolution for hadronic energy, so the E™ 1SS resolution is dominated by the 
latter. The events in both control samples are reweighted to reproduce the diphoton transverse 
energy distribution in the candidate 77 sample and, therefore, the transverse energy of the 
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hadronic recoil against the diphoton system; the reweighting factors range from 0.3 to 1.7. 
The shapes of the E™ 1SS distributions obtained from both control samples are identical within 
the statistical and systematic uncertainties and are used to predict the E™ 1SS background by 
normalizing to the number of events with E™ 1SS < 20 GeV in the candidate sample. 

The second background comes from events with real E™ 1SS . It is dominated by events with a 
genuine or fake photon and a W that decays into a neutrino and an electron, with the latter 
misidentified as a photon. This background is referred to as the electroweak background. An 
electron is misidentified as a photon if it satisfies all the photon selection criteria but has no 
matching hit pattern in the pixel detector. To model this background, an ey candidate sample 
is defined, selected the same way as the 77 sample but requiring at least one electron and at 
least one photon instead of at least two photons. This sample is thus enriched with W7 and W 
plus jets events (with real E™ 1SS ) similar to those in the candidate 77 sample. The probability, 
/ e ->7, to misidentify an electron as a photon is measured with Z — » ee events to be (1.4 ± 0.4)%. 
The e7 sample is then weighted by fe-^y/ (1 — fe^y) in order to estimate the contribution of 
this background to the signal sample. 

Finally, high-energy muons from cosmic rays or beam halo can deposit a large amount of en- 
ergy in the ECAL, leading to events with two photon candidates and E™ 1SS . These events are 
suppressed to a negligible level by the jet requirement described above. 

The E™ ss distribution in the 77 sample is shown in Fig. [lj together with the estimates of the 
electroweak background and the total background using the QCD prediction from Z — > ee and 
their uncertainties. Events with E™ 1SS < 20 GeV have negligible SUSY signal contribution (0.022 
± 0.009 event for the GGM SUSY sample point described below). The E™ 1SS distributions for 
the ff and ee samples are then scaled so that their integrals below 20 GeV match that of the 
77 sample minus the estimated electroweak contribution. Their integrals above 50 GeV give 
predictions of the QCD background. 

Table [l] summarizes the number of 77 events observed and the numbers of background events 
expected with Eij? iss > 50 GeV. The statistical and systematic uncertainties on the background 
predictions due to reweighting and normalization are shown separately. One event is observed, 
while the total background is expected to be 0.53 ± 0.37 (1.71 ± 0.64) events using the ff (ee) 
samples and including the electroweak background. These two consistent estimates are aver- 
aged, using log-normal distributions as probability density functions while taking into account 
the common component from the electroweak background and the correlated uncertainty due 
to the normalization, to obtain a prediction of 1.2 ± 0.4 background events. An additional 
conservative systematic uncertainty of 0.7 events is assigned by taking the largest difference 
between the average and the individual measurements, resulting in a total background uncer- 
tainty of 0.8 events. 

The efficiency for SUSY events to satisfy the selection criteria is determined by applying cor- 
rection factors derived from the data to the MC simulation of the signal. Since there exists no 
large clean sample of isolated photons in the data, we rely on similarities between the detector 
response to electrons and photons to extract the photon efficiency. The difference between the 
efficiencies for electrons and photons according to the MC simulation is 0.5%. The ratio of the 
electron efficiency from Z — > ee events that pass all photon identification criteria (except for the 
pixel match) to the corresponding electron MC efficiencies gives an MC efficiency scale factor 
of 0.967 ± 0.015. This scale factor is applied to the efficiency to identify a photon from NLSP 
decay obtained with the MC simulation. The error on the scale factor also includes possible 
systematic effects from multiple interactions per bunch crossing (pile-up), estimated to be less 
than 1%. The additional uncertainty on the efficiency of the requirement for no pixel match for 
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Figure 1: Eij? iss distribution for 77 data, including the jet requirement, compared with back- 
grounds and a possible GGM SUSY signal. The solid circles with error bars represent the data. 
The double-hatched blue band represents the contribution of the electroweak background. The 
single-hatched red band shows the sum of the electroweak background with the QCD E™ ss 
prediction obtained from the Z — > ee sample. The widths of the bands correspond to the sum 
of the statistical and systematic uncertainties on the backgrounds. The prediction of the GGM 
SUSY sample point described in the text is shown in the plot as the solid line histogram. 



Table 1: The number of events with E™ 1SS > 50 GeV from the 77 event sample as well as the 
predicted number of background events with E™ 1SS > 50 GeV using either the fake-fake events 
or the Z — > ee data. 



Type 


Number of 


Stat 


Reweight 


Normalization 




events 


error 


error 


error 


77 events 


1 








Electroweak background estimate 


0.04 ± 0.03 


±0.02 


±0.0 


±0.01 


QCD background estimate (ff ) 


0.49 ± 0.37 


±0.36 


±0.06 


±0.07 


QCD background estimate (ee) 


1.67 ±0.64 


±0.46 


±0.38 


±0.23 


Total background (using ff) 


0.53 ± 0.37 








Total background (using ee) 


1.71 ±0.64 








Combined total background 


1.2 ±0.8 








Expected from GGM sample point 


8.0 ±1.7 
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a photon is estimated by varying the amount of tracker material in the detector simulation and 
is found to be equal to 0.5%. Other sources of systematic uncertainties in signal yield include 
the uncertainty on the integrated luminosity (11%) fl29|, parton distribution function (PDF) un- 
certainty (10-40%), and renormalization scale uncertainty (10-20%), depending on the SUSY 
masses. 

As a cross check of the analysis method, in particular the determination of the QCD background 
contribution, the procedure is applied to a data sample selected in the same way as the 
sample described above, except that there is no requirement on the number of jets. The check 
consists of estimating the QCD background, events with no true E™ 1SS , in a sample populated 
with both QCD and events from the two dominant SM processes, W7 and W plus jets, which 
have true E™ 1SS . An excess in the observed number of events with large E™ 1SS over the estimated 
QCD background should be consistent with the yield from the W processes. The E™ 1SS spectrum 
of the e7 events is shown in Fig. [2] and exhibits a clear deviation from the predicted QCD 
background alone but agrees well with the sum of the QCD and W backgrounds. 
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Figure 2: The E™ ss distribution of the candidates (solid circles with error bars) is compared 
to the QCD expectation for this spectrum from Z — > ee events (solid histogram). The W7 
and W plus jet (with the jet misidentified as a photon) contributions are shown as dashed and 
dash-dotted histograms, respectively. The hatched red band represents the sum of all expected 
background components, and its width corresponds to the sum of the statistical and systematic 
uncertainties on the backgrounds. 



This measurement and the estimated acceptance times efficiency are used to set upper limits 
on the gluino and squark production cross sections, employing a Bayesian method described 
in Ref. [30 j. Both the log-normal and gamma priors are used to integrate over nuisance param- 
eters in order to incorporate uncertainties on the total background rate, integrated luminosity, 
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and total acceptance times efficiency. The observed 95% confidence level (CL) cross section lim- 
its, shown in Fig.[3j vary between 0.3 and 1.1 pb for a neutralino mass of 150 GeV as a function 
of squark and gluino masses. The variation in the limits is due to the dependence of the photon 
isolation efficiency on the number of jets in the event. 

To illustrate the limit-setting procedure, a particular signal point, with a squark mass of 720 GeV, 
gluino mass of 720 GeV, and neutralino mass of 150 GeV, is considered. The NLO signal cross 
section for this point is 1.04 pb, the efficiency times acceptance is 0.203 ± 0.004 (stat) ± 0.008 
(syst), and the PDF and scale uncertainties are each 13%. The upper limit obtained from this 
measurement is 0.585 pb, while the expected cross section upper limit at the 95% CL is 0.628 
pb. The expected number of events for this SUSY sample point for an integrated luminosity of 
36 pb -1 is 8.0 ± 1.7 events. For comparison, the total NLO SUSY production cross section at 
the Tevatron for the same model point is 0.3 fb, leading to an expectation of less than half of an 
event for analysis dJ. 

The benchmark GGM model used in this Letter can be used to interpret these cross section lim- 
its as lower limits on squark and gluino masses. For each mass point we compare the predicted 
cross section with the measured upper limit and claim exclusion if the former is greater than 
the latter. The exclusion contours for three different choices of neutralino mass are shown in 
Fig.|4j together with the expected exclusion limit for a 150 GeV neutralino mass. The predicted 
cross section has uncertainties due to the choice of PDF and scale. We varied the predicted 
cross section by one standard deviation of its uncertainty to ascertain the impact on the exclu- 
sion region. The results for 150 GeV neutralino mass are presented as a shaded band around 
the central value of the exclusion contour. 

In summary, a search for evidence of GGM SUSY production in events that contain two or more 
high transverse energy isolated photons, one or more jets, and large E™ ss is presented. No such 
evidence is observed, leading to upper limits on the GGM SUSY cross section between 0.3 and 
1.1 pb at the 95% CL across the parameter space of a benchmark model, setting the world's best 
direct lower limits on squark and gluino masses. 
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